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Many dynamic systems containing a large number of modes can benefit from adaptive 
control techniques, which are well suited to applications that have unknown parameters and 
poorly known operating conditions. In this paper, we focus on a model reference direct 
adaptive control approach that has been extended to handle adaptive rejection of persistent 
disturbances. We extend this adaptive control theory to accommodate problematic modal 
subsystems of a plant that inhibit the adaptive controller by causing the open -loop plant to 
be non-minimum phase. We will augment the adaptive controller using a Residual Mode 
Filter (RMF) to compensate for problematic modal subsystems, thereby allowing the system 
to satisfy the requirements for the adaptive controller to have guaranteed convergence and 
bounded gains. We apply these theoretical results to design an adaptive collective pitch 
controller for a high-fidelity simulation of a utility-scale, variable-speed wind turbine that 
has minimum phase zeros. 


I. Introduction 

M any plants, including wind turbines and rotor craft, have modal subsystems that interfere with the ability of an 
adaptive controller to perform in a stable manner. In this paper, we extend our model reference adaptive 
control theory [l]-[4] to accommodate modal subsystems of a plant that inhibit the adaptive controller. In particular, 
residual modes that have non-minimum phase zeros can inhibit the controller by interfereing with the almost strict 
positive real condition. This theory applies to linear time invariant systems that can be diagonalized or placed into 
modal form. Linear flexible structures are good examples of appropriate applications of this theory. 

The controller modification will use the idea of Residual Mode Filters (RMF) introduced for fixed gain controllers 
in [5], In this paper, the RMF will be used to eliminate the effect of modes that prevent the almost strict positive 
realness of the overall system. These modes have non-minimum phase zeros. 

II. Direct Model Reference Adaptive Control with Rejection of Persistent Disturbances 

In this paper, we focus on the direct adaptive control (DAC) approach developed in [l]-[2] and extended in [3]-[4] 
to handle adaptive rejection of persistent disturbances. The plant used in this section of the paper will be modeled by 
the linear, time -invariant, finite -dimensional system: 

[x p =Ax p +Bu p 
\y P =Cx p ;x p (0) 


+ Tm, 


= X n 


( 1 ) 
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where the plant state, X p (t) , is an N p -dimensional vector, the control input vector, u p (t) , is M-dimensional, and 
the sensor output vector, y p (t), is P-dimensional. The disturbance input vector, u D (t) , is M D -dimensional and 
will be thought to come from the Disturbance Generator: 


[u D =&Z D 

\z D = F z D ;z D (0) = z 0 


( 2 ) 


where the disturbance state, z D (t) , is N D -dimensional. All matrices in Eqs. (l)-(2) have the appropriate compatible 
dimensions. Such descriptions of persistent disturbances were first used in [6] to describe signals of known form but 
unknown amplitude. Equation (2) can be rewritten as in [3] in a form that is not a dynamical system, which is 
sometimes easier to use: 


j u D =® z d 

{ Z D = Ffio 


(3) 


where (j) n is a vector composed of the known basis functions for the solution of u D — &z D , i.e., (j) D are the basis 
functions which make up the known form of the disturbance, and Lisa matrix of dimension ND by dim (</> D ) . The 
method for rejecting persistent disturbances used in this paper requires only the knowledge of the form of the 
disturbance, the amplitude of the disturbance does not need to be known, i.e. ( L ,& ) can be unknown. In this paper, 
we will rejecting step disturbances of unknown amplitude which can be represented in the form of Eq. (3) as 
</> D = 1, with ( L,Q) unknown. 

In [7], as with much of the control literature, it is assumed that the plant and disturbance generator parameter 
matrices, , are known. This knowledge of the plant and its disturbance generator allows the 

Separation Principle of Linear Control Theory to be invoked to arrive at a State -Estimator based, linear controller 
that can suppress the persistent disturbances via feedback. In this paper, we will not assume that the plant and 
disturbance generator parameter matrices (A,B,C,T,@) are known. Instead, we assume that the disturbance 
generator parameter F from Eq. (2) is known, i.e., the form of the disturbance functions is known. In many cases, 
knowledge of F is not a severe restriction, since the disturbance function is often of known form but unknown 
amplitude. 

Our control objective will be to cause the plant output y p (t) to asymptotically track the output of a known 
reference model y m ( t) . The reference model is given by: 


{ x = A x + B u : x (0) = x 

m m m m m ’ m V s 

y m ^ m ^ m 


m 

0 


(4) 


where the reference model state, X m (t) , is an N m -dimensional vector. The reference model output, y m ( /) , must 
have the same dimension as the plant output, y p (t) . The excitation of the reference model is accomplished via the 
vector, u (t) , which is generated by: 


u = F u ; u (0) = (5) 

m m m ’ m V ' 0 V' 

It is assumed that the reference model is stable and the model parameters ( A m ,B m ,C m ,F m ) are known. 

The desired control objective is for the output of the plant to asymptotically track the output of the reference 
model. We define the output error vector as: 


e y =y p -y m 


(6) 
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where h m ,h u ,h e , and h D are arbitrary, positive definite matrices. The adaptive controller is specified by Eq. (8) 
with the adaptive gain laws given by Eq. (9). See [4] for an analysis of the stability of this controller and the proof of 
the following result: 

Theorem 1: Suppose the following are true: 

(1 ) All u m are bounded (i.e., all eigenvalues of F m are in the closed left-half plane); 

(2) The reference model (4) is stable (i.e., all eigenvalues ofA m are in the open left-half plane); 

(3) (j) D is bounded (i.e., all eigenvalues of F are in the closed left-half plane and any eigenvalues on the 
imaginary axis are simple); 

(4) (A, B, C) is Almost Strict Positive Real (ASPR), i.e., CB> 0 and the open-loop transfer function is 
minimum phase. 

Then the adaptive gains G m , G u , G e , and G D , remain bounded and asymptotic tracking occurs, i.e., 

e =v —V >0. 

y ' p oo 


III. Residual Mode Filter Augmentation of Adaptive Controller 

In some cases the plant in Eq. (1) does not satisfy the ASPR requirement of the adaptive controller. Instead, there 
may be be a modal subsystem that inhibits this property. In [8]-[9] we developed theory to augment the adaptive 
controller using a Residual Mode Filter (RMF) to compensate for the troublesome modal subsystem, or the Q 
modes. Here we present new theory for augmented adaptive model reference controllers using RMF. 

Let us assume that Eq. (1) can be partitioned into the following modal form: 
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Define X p = 


X 

, A = 

A 

0 " 

, B n = 

~ B~ 

, r n = 

~ r " 

, C T = 

~ c~ 

_X e _ 

’ P 

0 

Aq_ 

’ p 

B a_ 

’ p 

_ £ T q_ 

’ p 

c a_ 


generator 


I z D = Fz d 
[ u D = 0z D 


and the disturbance 


The output tracking error and control objective remain as in Eqs. (7)-(8), i.e.. 


e y =y p —y m — y-y — >0. However, now we will only assume that the subsystem ( A,B,C ) is Almost Strictly 

Positive Real, rather than the full un-partitioned plant (A p , B p , C p and the modal subsystem (Aq,Bq, C 0 ) will be 

known and open-loop stable, i.e., A 0 is stable. Also note that this subsystem is directly affected by the disturbance 

input. In [10] it was shown that a system ( A,B , C ) is ASPR when CB is positive definite and the open-loop system 

P(s) = C(sI — A ) 1 B is minimum phase. So, in summary, the actual plant has an ASPR subsystem and a known 
modal subsystem that is stable but inhibits the property of ASPR for the full plant. Hence, this modal subsystem 
must be compensated or filtered away. 

We define the Residual Mode Filter (RMF): 


l*e =A Q x Q +B Q u p 


y 0 = c oF 


(id 


Q~Q 


And the compensated tracking error: 


e y = e y~yQ 


( 12 ) 


Now let e ( j = X Q — X Q and use substitution to obtain the error system: 


—AqCq sYqU d 


(13) 


Consequently we can use Eq. (11) and the definition of e () to write Eq. (12) as: 

e y =(y P -y m )-c.i Q 

~ (C',x + CqXq — Cx m (FqXq + Cq€q J (14) 

= CAx - C Q e Q 

As in [l]-[2], we define the ideal trajectories X„ and u* as a linear combination of the plant states, the control 
inputs, and the disturbance inputs. Note that this is done only for the ASPR subsystem (A,B, C j as follows: 

fx, = Ax* +Bu t +Yu D \x t ( 0 ) = x 0 

^ (15) 
[T* ^ Cx t = y m 


with J X * _ ^1 Pm + S\2 U m + ^13 Z £> 

[w* = S 2 ]X m +S 22 u m + S 2 - i z D 

Note that the ideal output y„ matches the reference model output y . If such ideal trajectories exist, they will 
produce exact output tracking. The ideal trajectories can be made to match the reference model in Eqs. (4)-(5) with 
the following model matching conditions'. 

4 

American Institute of Aeronautics and Astronautics 



( 16 ) 


} S*F = AS* + BS 2 +to 

W = C B 


where 5, = [S u S l2 

s; = [s;, s : 2 s:,] 

F = 

0 

1 

F ,n 

0 0 

1 

, r = 

1 

0 O 

1 

C m = [C m 0 0], 



l_o 

0 

F 


L r J 


The model matching conditions given in Eq. (16) are necessary and sufficient conditions for the existence of 
ideal trajectories. Solutions to these matching conditions must exist for later analysis, but explicit solutions need 
never be known for the adaptive controller design. The model matching conditions are know to be uniquely solvable 
when CB is nonsingular [10]. 


Let] 


Ax = x-Xt, 

Am = u p - m* , then from Eqs. (14) and (15) we can write: 
Ay = e v 


At = A Ax + BAu 
Ay = CAx-C Q e Q 


(17) 


This system can be rewritten as: 


Ay 


0 

Ay 


~B~ 


0 


Ay 


= 



+ 


A u + 


u n = A 


}q\ 


0 A, 

_ e a\ 


0 


1 

l 

c*} 

n 

(O 

1 

LJ 

_ e Q_ 


+ Z? A u + sY' qU q 


r 1 

Ay 


Ay 

a y = [c -c Q \ 

_ e Q\ 

= c 

/G_ 


(18) 


Now we have the following: 



r_ 

"A 0 " 


~B~ 

A 

— r l 

Lemma: 

A = 


M = 


,C =\C -C 0 \ 


V 

1 

O 

to 

1 


0 

’ L Q J 

J 


Proof: 


: cs = [c -cj 


= CB > 0 and 


P{s) = C{sI -A) { B 


= [c -c Q ] 


(si -Ay' 

0 


0 

(. si-A 0 y x 


= C(sl - A) 1 B = P(s) 

is minimum phase. End of proof. 

So there exists G e such that (A c = A + B G e C ,B ,C) is Strictly Positive Real (SPR) when ( A,B,C) is ASPR. 
Consequently, as is well known from the Kalman- Yacubovic Theorem, there exists P ,Q> 0 such that 
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( 19 ) 


\a t c p+pa c =-q 
\pb = c t 


We now write the augmented adaptive control law using RMF: 


Up =G m x m +G u u m +G e e y +G D </> L 

e y s -v, - y q 


x Q = AqXq + B Q u p 


y Q = Gqx c 


with adaptive gains given by: 


G m =~e y x m h m ;h m >0 
G„ =-e y u m h H ;h u >0 
G e = —e y e'h e \h e >0 
Go ~ oh o -ho ^ 0 


Finally, we have the following stability result: 

Theorem 3: In Eq. (10), let (A, B, C ) be ASPR, A (J be stable, and </) r> be bounded. Then the augmented adaptive 
controller using RMF in Eqs. (20)-(21) produces e = y p — V m and e Q ultimately bounded into a ball of radius 

(l + V^inax) 

R* = £ . — M v with exponential rate and bounded adaptive gains (G m ,G U ,G e ,G D ) . 

U \ P min 

Proof: From Eq. (20), we have U = G m X m + G u U m + G e e y + G D <j> D , so we can write 


R t =£ 


A u = u p - u* 

= \G x + G u + 

L m m u m 

= GX+AGj 1 
A G u =G U -Si, 
A G m - G m - Si 
A G e =G e -G* e 


G e e y +G D (f > D ] [A 21 T m + S 22 u m + S 2i L<j> D ] 


AG d =G d -(S* 23 L) 

where j AG = G-G„= [AG m A G u A G e A G d \ 
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Then 


£, = A £ + BAu = A c £ + Bw + sT Q u D 
e v =C^ 


( 22 ) 


with C, = 


Ax 


, w = AGr/,A c =A + BG e C . 


From Eq. (21), we can see that 


G = AG = -e n T h;h = 


K o o 

o h u o 

o o h e 

0 


0 

0 

0 


0 


0 hr 


>0 


(23) 


Since ( A,B,C) is ASPR, and by the lemma, so is (A.B,C), we can use the following result from [11] where 
v = T () u d is bounded because the disturbance u D —L(f> D is bounded. 

Result: Consider the nonlinear, coupled system of differential equations, 

g = A c g + B (g(0 - G * ) r/+ sv 

e y =Cg (24) 

G(t) = -e v r/ T h - aG(t) 


where O is any constant matrix and h is any positive definite constant matrix, each of appropriate dimension. 
Assume the following: 
i) the triple (A, B ,C) is SPR, 


ii) there exists M K > 0 such that 


< M k , using the trace norm, 


iii) there exists M v > 0 such that S U p| | V'C — M v 


t> o 


iv) there exists a > 0 such that a < 


qm m and 


2 Pn 


v) h satisfies \\h 1 < 


^2 


sM v 

\ClM K j 


minimum eigenvalue of Q in the system ^ 

pb=c t 


, where p mi „, p max are the minimum and maximum eigenvalues of P and q min is the 

\a t c p+pa c =-q 


R*=£ 


( 1 + A 

1 P max ) 

a sl 

P min 


M ,, with £>0. 
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From this result, we have ^ is ultimately bounded into the ball of radius R>, which leads to 
e = y —y m — y p — y* =CA and e 0 ultimately bounded as well. Therefore G = G* + AG is bounded, as desired. 
End of proof. 

Consequently, the radius of the error ball R, =£— . 7 M v is determined by the size of s, which is 

WA nin 

related to the amount of disturbance leakage into the Q modes. It can be seen that, when there is no leakage of the 
disturbance into the Q modes (s = 0), the convergence is asymptotic to zero. 

IV. Design of Wind Turbine Controller 


A. Wind turbine simulation 

Rated wind speed is the velocity at which maximum power output, or rated power, of a wind turbine is achieved. 
If a turbine is allowed to operate in an uncontrolled manner, in conditions where the wind speed is above the rated 
wind speed, the power output would increase in proportion to the cube of the wind speed, resulting in overheating of 
the generator and the power electronics system. Additionally, high wind speeds result in larger aerodynamic forces 
on the machine, causing fatigue and possible system failure. Hence, power output of a turbine must be held constant 
for wind speeds at and above the turbine’s rated wind speed. This wind turbine operation area is called Region 3 
[12], For variable-speed turbines operating in Region 3, a constant torque is applied at the generator, and the turbine 
rotational speed is maintained at the desired value through the use of collective blade pitch. 

Wind turbine control problems can benefit from adaptive control, which are well suited to nonlinear applications 
that have unknown modeling parameters and poorly known operating conditions. The main nonlinearities in a wind 
turbine model come from the nonlinear aerodynamic loads on the turbine. Creating an accurate model of the 
dynamic characteristics of a wind turbine is expensive and extremely difficult, if not impossible. Additionally, wind 
turbines operate in highly turbulent and unpredictable conditions. These complex aspects of wind turbines make 
them attractive candidates for the application of adaptive control methods. 

The Controls Advanced Research Turbine (CART2) is a two-bladed, upwind, active-yaw, variable-speed 
horizontal-axis wind turbine located at the National Wind Technology Center (NWTC) in Golden, Colorado. This 
machine is used as a test bed to study aspects of wind turbine control technology for medium-scale machines [1], 
The pitch system on the CART uses electromechanical servos that can pitch the blades up to ±18 deg/s. In Region 3, 
the CART uses a conventional variable-speed approach to maintain rated electrical power, which is 600 kW at a 
low-speed shaft [LSS] speed of 41.7 RPM and a high-speed shaft [HSS] speed of 1800 RPM. Power electronics are 
used to command constant torque from the generator and full-span blade pitch controls the turbine speed. 

The CART has been modeled with the Fatigue, Aerodynamics, Structures, and Turbulence Codes (FAST) as a 
combination of rigid and flexible bodies connected by several degrees of freedom (DOFs). The DOFs can be turned 
on or off individually for analysis purposes by setting a switch in the input data file. The FAST Code is a 
comprehensive aeroelastic simulator capable of predicting both the extreme loads and the fatigue loads of two- and 
three-bladed horizontal axis wind turbines [13]. FAST uses Kane’s method to set up equations of motion that are 
solved by numerical integration. The aerodynamic forces and moments along the turbine blade are calculated in 
FAST by the AeroDyn subroutine package [14]. 

B. Region 3 Collective Pitch Controller 

In this section, we use the theory from section II to design region 3 controllers for the FAST simulation of the 
CART2. The parametric information for the FAST simulator as we configured it is available from [15], For the 
adaptive controller with a reference model, the control objective is to track the model, to reject wind disturbances, 
and to operate in the presence of flexible modes using collective blade pitch. The model is designed such that output 
tracking results in plant generator speed regulation at 1800 rpm. An adaptive controller without a reference model as 
described in ref. [9] has the objective of regulating the generator speed instead of tracking a model. The inputs to the 
FAST plant are generator torque, blade pitch angle, and nacelle yaw. The FAST simulator can be configured to 
output many different states or measurements of the plant, such as generator speed and low speed shaft velocity. In 
this study, the yaw is assumed fixed, so that the major component of the wind inflow is normal to the rotor. In 
addition, the generator torque is assumed constant in region 3. Collective blade pitch is the only controller output. 
Turbine rotational speed, measured on the low-speed shaft side of the gearbox, is the only plant output used by the 
region 3 controllers. A classical PI collective pitch controller (the baseline PI pitch controller) has been implemented 
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and tested in the FAST simulator and a similar controller was field tested on the CART2 [16]-[17] for validation of 
the codes. We use the baseline PI pitch controller as a basis for comparison with the adaptive pitch controller. 

The adaptive collective pitch controllers designed for this paper replace the baseline PI pitch controller in the 
FAST simulator. The uniform wind disturbance, without shear, across the rotor disk of a turbine can be modeled as a 
step disturbance [18]. Hence, to improve controller performance and reduce loads due to changes in wind speed, we 
design the adaptive collective pitch controller to reject step disturbances of unknown amplitude. The control 
objectives are accomplished by collective blade pitch. 

A control law of the form given in Eq. (21) with gains specified by Eq. (20) is used to design the model 
reference adaptive collective pitch controller. A step function is used as the disturbance generator function, i.e., 
(f> D = 1 from Eq. (2). Recall that the amplitude of the disturbance function does not need to be known. This 

adaptive controller is implemented in Simulink™ for the FAST simulation of the CART2. The reference model 
chosen was a linearized model including the generator state and the tower state. The poles corresponding to the 
tower first fore-aft bending mode were moved from -0.0722±5.4999/ to -0.2722±5.4999/ using output feedback. The 
adaptive controllers’ gains were tuned to track the reference model, minimize the generator speed error, and reject 
the wind disturbances. The gains used in the model reference adaptive controller are: h u = 0.01, h m = 0.1, h e = 8.5 and 
h D = 0.3. The gains used in the regulating adaptive controller are:, h e = 6.5 and h D = 0.3. 

In a previous study 8 , the blade first flap-wise bending DOF when enabled along with the drive-train mode, were 
shown to cause the adaptively controlled turbine to saturate the blade pitch actuators after 60 seconds of steady 
uniform wind inflow at 18 mps, causing the turbine to enter over-speed condition. The baseline PI controller 
remained within acceptable operating limits, but the generator speed tracking was unacceptable at wind speeds 
above 17 mps. 

A stability condition of the adaptive controller is an ASPR plant, i.e., CB positive definite and the open-loop 
transfer function of the plant has no non-minimum phase zeros. When the turbine simulation is trimmed at a wind 
speed of 18 mps with the generator and drive-train DOFs enabled to obtain a linear model of the plant, the open-loop 
transfer function of the linearized plant model has two non-minimum phase zeros at 0.01 1 1+5.499/. Hence the plant 
does not satisfy the ASPR condition for the adaptive controller. We use the theory from section II above to design a 
Residual Mode Filter to remove the non-minimum phase modes from the plant to restore the ASPR property of the 
plant. 

A Residual Mode Filter is designed from the linear model by first converting the linear system to a modal 
system. We partition the modal system into two subsystems, one minimum phase subsystem and a second stable 
subsystem with two non-minimum phase zeros. The second subsystem contains the Q-modes, so it is used as the 
Residual Mode Filter given in Eq. (1 1) to augment the adaptive controller to remove the plant’s non-minimum phase 
modes. The transfer function for the RMF is 

, 790.52s -38.32 

T(S ) = — . ( 22 ) 

s 2 + 0.02s + 430.30 

The RMF is placed in a loop around the controller in the Simulink™ model of the turbine. The controller output 
is fed to both the plant and the RMF. The RMF output is subtracted from the plant output before it is passed to the 
controller. The RMF removes the modes from the plant output that inhibit the ASPR property. Previous results from 
[14] showed that fixed gain controllers could be augmented with RMFs, so the same RMF is added to the original 
baseline PI pitch controller. An advantage of RMF augmentation is that it requires no modification to the control 
laws or control gains. The RMF is designed to remove only the Q-modes, whereas a low-pass filter attenuates or 
removes frequencies above a certain value. Next we present simulation results using these controllers. 

C. Simulation Results 

The baseline and model reference adaptive controllers augmented using the RMF with transfer function given by 
Eq. (22) were compared in simulation. Simulations were run with turbulent wind inflow to excite the nonlinear 
turbine model. The turbulent wind primarily results in region 3 operation, see fig. 1 . All of the degrees of freedom 
for the simulation were enabled. The resulting generator speed from the simulations are shown in fig. 2. The PI 
controller exceeds the generator speed limit at 30 seconds. Figure 3 shows the comparison of the model reference 
adaptive controller with the adaptive controller with no model (the regulating adaptive controller). Both adaptive 
controllers have comparable results. Preliminary studies of the tower bending modes did not show significant 
changes when the model was used with the adaptive controller. 
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Time (s) 

Fig. 1. Turbulent wind inflow used in simulations. 



Time (s) 

Fig. 2. Baseline PI with RMF and Reference Model Adaptive RMF with turbulent 
wind inflow. 
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Time (s) 

Fig. 3. Generator speed for RMF adaptive controllers with and without a reference model. 


V. Conclusion 

We have proposed an augmented model reference adaptive controller with adaptive rejection of persistent 
disturbances using a residual mode filter. The RMF is used to accommodate problematic modes in the system that 
inhibit the adaptive controller, in particular the ASPR condition. This new theory accounts for adaptive model 
tracking. This theory is used to design a controller for a utility scale wind turbine that outperformed the baseline PI 
controller. Future work will include investigating different turbine reference models to determine if they can have a 
positive affect on the loads experienced by the turbine. 
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